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We investigate the properties of light (anti)nuclei and (anti)hypertriton production in Pb-Pb 
collisions at ^snn = 2.76 TeV, based on the parton and hadron cascade and dynamically constrained 
phase-space coalescence (PACIAE-l-DCPC) model. We found that the yields of light (anti)nuclei 
and (anti)hypertriton strongly depend on the centrality, i.e., their yields decrease rapidly with the 
increase of centrality bins; but their yield ratios are independent of centrality. The results of our 
theoretical model are well consistent with ALICE data. Furthermore, we found that the integrated 
yields of (anti)nuclei per participant nucleon increase from peripheral to central collisions and more 
rapidly with increasing mass number. The transverse momentum distributions of the d, d, 

\H, ^He and ^He are also discussed in the 0-10% most central Pb-Pb collisions. The coalescence 
parameters Ba of light (anti)nuclei and (anti)hypernuclei are analyzed. 


I. INTRODUCTION 

Antimatter production has received considerable at¬ 
tention in particle and nuclear physics, astrophysics, cos¬ 
mology and other fields of modern physics, since Dirac 
predicted the existence of negative energy states(i.e., an¬ 
timatter) of electrons in 1928. The antiprotons [l| were 
discovered in 1955, and antineutrons were discov¬ 
ered in 1956, followed by a series of antimatter nuclei 
including anti-deuterons, anti-triton, the anti-helium-3 
in scientific experiments dl. The hot and dense mat¬ 
ter forms in ultra-relativistic heavy-ion collisions, which 
contains roughly equal numbers of quarks and anti¬ 
quarks 0, is similar to the fireball environment in the 
initial stages of the Big Bang.Thus, the ultra-relativistic 
eavy-ion collisions provide a good experimental condi¬ 
tion Q to study the production of light (anti)nuclei and 
the evolution of the early universe. 

In 1995, the first production of anti-hydrogen atoms 
was detected Q, and then which was trapped success¬ 
fully with a confinement time of 172 ms by the AL¬ 
PHA at the European Organization for Nuclear 
Research (CERN). Moreover, the ALICE Collabora¬ 
tion has also published its preliminary d yield of ^ 
6 X 10“® measured in the proton-proton collisions at 
yi = 7 TeV m- At BNL(the Brookhaven National 
Laboratory), the STAR Collaboration has reported their 
measurements of TiL [l^ and ^He Qin Au-Au collisions 
at the top RHIC energy, respectively. 

In theory, firstly the nucleons and hyperons are usu¬ 
ally calculated with some selected models, such as the 
transport model. Then, the light nuclei (anti-nuclei) 
are studied by the reasonable hadron final-state coales¬ 
cence models jjMa, such as the phase-space coales¬ 
cence model [imj and/or the statistical model [^[^. 
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etc. For example, the production of light nuclei (hyper- 
nuclei) in the Au-Au and Pb-Pb collisions at relativistic 
energies was described theoretically by the coalescence 
-I- blast-wave method 23l| an d the UrQMD-hydro hybrid 
model + thermal model [24| , respectively. 

Recently, we have proposed an approach studying 
the light nuclei (anti-nuclei) production in the relativis¬ 
tic heave ion collisions by the dynamically constrained 
phase-space coalescence model(DCPC) [^, which is 
based on the final state hadronic generated by a par- 
ton and hadron cascade model PACIAE [2^. Using 
this method, the light nuclei (anti-nuclei) yields, trans¬ 
verse momentum distribution, and rapidity distribu¬ 
tion in non-single diffractive proton-proton collisions 
at ^/s = 7 TeV [2^ are predicted. The light nuclei 
(anti-nuclei) and hypernuclei (anti-hypernuclei) produc¬ 
tions , their centrality dependence [2^ and their mass 
number scaling property (29| in the Au-Au collisions at 
^snn = 200 GeV are also investigated. Moreover,The 
energy dependence of the ratio for antiparticle to par¬ 
ticle is studied in the high energy proton-proton colli¬ 
sions [ 3 ^. In this paper, we will use this method to 
investigate the light nuclei (anti-nuclei) and hypernuclei 
(anti-hypernuclei) productions and the properties in the 
Pb-Pb collisions at .^/snn = 2.76 TeV. 

The paper is organized as follows: In Sec. II, we briefly 
introduce the PACIAE and DCPC model. In Sec. Ill, our 
calculated results of light nuclei (anti-nuclei) and hy¬ 
pertriton (anti-hypertriton) are presented, such as their 
yields, ratios, the transverse momentum distributions, 
and are compared with the ALICE data. In Sec. IV, a 
short summary is given. 


II. MODELS 

The PACIAE model [ 2 ^ is based on PYTHIA6.4 [3l| 
and is designed mainly for nucleus-nucleus collisions. In 
the PACIAE model, the process is decomposed into four 
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steps. Firstly, the nucleus-nucleus collision is decom¬ 
posed into the nucleon-nucleon (A^A^) collisions according 
to the collision geometry and NN total cross section. 
Each NN collision is described by the PYTHIA model 
with the string fragmentation switches off and the di¬ 
quarks (antidiquarks) randomly breaks into quarks(anti- 
quarks). So the consequence of a NN collision is a 
partonic initial state composed of quarks, anti-quarks, 
and gluons. Provided all NN collisions are exhausted, 
one obtains a partonic initial state for a nucleus-nucleus 
collision. This partonic initial state is regarded as 
the quark-gluon matter(QGM) formed in relativistic 
nucleus-nucleus collisions. Secondly, the parton rescat¬ 
tering proceeds. The rescattering among partons in 
QGM is randomly considered W the 2—>• 2 LO-pQCD 
parton-parton cross sections [32j| . In addition, a K 
factor is introduced here to account for higher order 
and non-perturbative corrections. Thirdly, hadroniza- 
tion happens after parton rescattering. The partonic 
matter can be hadronized by the Lund string fragmen¬ 
tation regim e 1311 and/or the phenomenological coales¬ 
cence model [26l| . Finally, the hadronic matter continues 
rescattering until the hadronic freeze-out (the exhaustion 
of the hadron-hadron collision pairs). We refer to [2^ for 
the details. 

With the final state particles provided by the PACIAE 
model, we can then calculate the production of light nu¬ 
clei (anti-nuclei)with the DCPC model. In quantum sta¬ 
tistical mechanics [s^, one cannot precisely define both 
position q = (x,y,z) and momentum p = (j)x,Py,Pz) of 
a particle in the six-dimension phase space because of 
the uncertainty principle AqAp^ h?. We can only say 
that this particle lies somewhere within a six-dimension 
quantum box or state with a volume of AqAp. A particle 
state occupies a volume of in the six-dimension phase 
space [ 13 . Therefore, one can estimate the yield of a 
single particle by defining an integral Yi = 
where H and E are the Hamiltonian and energy of the 
particle, respectively. Similarly, the yield of the N par¬ 
ticle cluster can be estimated as the following integral: 

dqidpi...dqNdpN . . 

h^N ■ 

In addition, equation o must satisfy the following 
constraint conditions: 



Too ^ ruinv ^ Too + Am, 
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^ ^ j]i,j = 1,2,...,N). (3) 


where 


N s 2 / Y X 21 

i—l ' ^ i—1 
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and Ei, pi(i=I,2,... ,N) are the energies and momenta of 
particles, respectively, mo and Dq stand for,respectively. 


the rest mass and diameter of light (anti)nuclei. Am 
refers to the allowed mass uncertainty, and qij = \qi — qj\ 
is the vector distance between particles i and j. Be¬ 
cause the hadron position and momentum distributions 
from transport model simulations are discrete, the inte¬ 
gral over continuous distributions in equation (HD should 
be replaced by the sum over discrete distributions. 


III. RESULTS AND DISCUSSION 

First we produce the final state particles using the 
PACIAE model [ 2 ^. In the PYTHIA simulations, we as¬ 
sume that hyperons heavier than A decayed already. The 
model parameters are fixed on the default values given 
in PYTHIA [3l[. However, the K factor as well as the 
parameters parj(l), parj(2), and parj(3), which are rel¬ 
evant to the hadrons production in PYTHIA, are given 
by fitting the ALICE data of p, p, A in Pb-Pb collisions 
at ^snn = 2.76 [H, . Specific details of this method 

is similar to the paper |28l |. The fitted parameters of 
K= 3(default value is 1 or 1.5 [3l|), parj(l) = 0.15(0.1), 
parj(2) = 0.38(0.3), and parj(3) = 0.65(0.4) are used to 
generate 1.833 x 10® minimum-bias events by the PA¬ 
CIAE model for the 0-20% centrality Pb-Pb collisions at 
^snn = 2.76 TeV with \y\ < 0.5 acceptances, as shown 
in Table ID 


TABLE I: The integrated yield dN/dy of particles at midra¬ 
pidity {\y\ < 0.5), for p,p and A in Pb-Pb collisions of 
x/snn = 2.76 TeV with 0-20% centrality. 


Particle type 

ALICE“ 

PACIAE 

P 

25.9 ± 1.6 

25.9 

P 

26.0 ± 1.8 

24.8 

A 

19.3 ± 1.4 

19.18 


“ The ALICE data are taken from Ref. [34,^^. 


Then, the integrated yields dN/dy of light (anti)nuclei 
d (d), ^He i^He), as well as \H {^H) are calculated by 
the DCPC model for each centrality bin of 0-5%, 5-10%, 
10-20%, 20-30%, 30-50% as well as 50-80%, as shown in 
Tab. El One can see from Table H that the yields dN/dy 
of light (anti)nuclei and (anti)hypertritons calculated by 
the DCPC model decrease (or increase) with the increase 
of centrality(or A^part); the yields of (anti)nuclei decrease 
with the increase of mass; and the yields of anti-nuclei 
are less than that of its corresponding nuclei. 

In order to facilitate comparison with the experimental 
data, the cumulate yield Yc is described as 





( 5 ) 


Where c is the value of centrality bins. The Fig. E effec¬ 
tively shows the cumulative yields Y/ of d, d, \H, AH, 
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TABLE II: The integrated yields dN/dy of d, d, ^He, ^He, \H and Ajf calculated by PACIAE+DCPC model in midrapidity 
Pb-Pb collisions of ^snn = 2.76 TeV. (Apart) is shown for each centrality bin. 


Centrality 

0% - 5% 

5% - 10% 

10% - 20% 

20% - 30% 

30% - 50% 

50% - 80% 

{^part) 

379 

323 

256 

186 

110 

36 


0.112 

0.0917 

0.0703 

0.0461 

0.0199 

0.00345 

T 

0.101 

0.0838 

0.0631 

0.0412 

0.0182 

0.00311 


3.84E-04 

2.95E-04 

2.13E-04 

1.12E-04 

4.99E-05 

5.13E-06 

HTl'’ 

3.29E-04 

2.51E-04 

1.82E-04 

9.44E-05 

4.30E-05 

4.29E-06 

IH^ 

4.47E-05 

3.40E-05 

2.33E-05 

1.16E-05 

5.94E-06 

5.67E-07 

W 

3.24E-05 

2.45E-05 

1.71E-05 

8.26E-06 

4.38E-06 

4.03E-07 


“ calculated with Am = 0.00042 GeV for d, d. 
calculated with Am = 0.00091 GeV for ^He, ^He. 
calculated with Am = 0.00040 GeV for \H, AA. 
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FIG. 1: The cumulative yields for light (anti)nuclei and 
(anti)hyper-triton as a function of centrality in midrapidity 
Pb-Pb collisions at ^snn = 2.76 TeV. The solid symbols 
are the experimental data points from ALICE [^, . The 

open symbols represent the outcome for our PACIAE-l-DCPC 
model. 


^He and ^He in different centrality bins Pb-Pb collisions 
at -y/sNN = 2.76 TeV. One can see in Fig. [U as well as 
Tab. m the yields of light (anti)nuclei and (anti)hyper- 
triton all decrease rapidly with the increase of centrality, 
presenting a fall as the index distribution. Meanwhile, 
the PACIAE+DCPC model results (the open symbols) are 
consistent with the ALICE experimental data 
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FIG. 2: Yield ratios of light anti-nuclei and anti-hypertriton 
(d, ^He, and ^H) to light nuclei and hyper-triton (d, 
^He, and \H) in midrapidity Pb-Pb collisions at = 

2.76 TeV, plotted as a function of centrality. Open symbols 
represent our PACIAE+DCPC model results. Solid symbols 
are the experimental data points from ALICE [3^ . IT7| , which 
statistical uncertainties are represented by bars and system¬ 
atic uncertainties are represented by open boxes. 

In Fig. [2l the yield ratios of light anti-nuclei and anti¬ 
hypertriton (d, ^He, and AA) to light nuclei and hyper¬ 
triton (d, ^He and \H), are given in different central¬ 
ity Pb-Pb collisions at = 2.76 TeV. To facilitate 

comparison, experimental results from ALICE [s^ 
are also given with the solid points. One can see from 
Fig. m the yield ratios of light anti-nuclei to light nu¬ 
clei and anti-hypertriton to hyper-triton from central to 
peripheral collisions remain unchanged, and their cor- 
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FIG. 3: The results (open points) of the transverse momentum distributions of light anti(nuclei) for PACIAE+DCPC model in 
the 0-10% most central Pb-Pb collisions at ^snn = 2.76 TeV with \y\ < 0.5 acceptances, calculated for (a) d and d, (b) 
and \H, and (c) ^He and ^He. 


responding values are respectively about 0.91, 0.85,0.72, 
although their yields decrease rapidly with the centrality 
as shown in Tab.lllland Fig.[T] The results obtained from 
our model are also in agreement with the experimental 
data from ALICE [s^, within error ranges. 
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FIG. 4: The integrated yield dN/dy as a function of 

A^part. The yields at midrapidity for p,p. A, A,d, d, \H, 
He,^ He are divided by Apart, normalized to the periph¬ 
eral collisions (50%-80%). The results are calculated by PA¬ 
CIAE+DCPC model in Pb-Pb collisions at Asnn = 2.76 TeV. 


In Fig.[3l the transverse momentum distributions of d, 
d, \H, ^He and '^He calculated were shown in the 


0-10% central Pb-Pb collisions at ^snn = 2.76 TeV with 
\y\ < 0.5 acceptances. Figures 3(a), 3(b) and 3(c) are the 
calculated the transverse momentum pt distributions for 
did), \Hi!LH), as well as ^He(^He), respectively. The 
pattern of light (anti)nuclei transverse momentum dis¬ 
tributions is similar to the ALICE 0113 of transverse 
momentum. In addition, you can see that the the peak 
of transverse distribution appears at similar range for 
nuclei and anti-nuclei situations. 

TABLE III: The results of the light (anti)nuclei average 
transverse momentum {pt) calculated by PACIAE+DCPC 
mode in the Pb-Pb collisions at ^snn = 2.76 TeV, compar¬ 
ing with the Au-Au collisions at y^SNN = 200 GeV [T7|. 


Particle type 

d 

d 

^He 


IH 


Pb-Pb(2.76 TeV) 

1.00 

1.05 

1.16 

1.24 

1.17 

1.22 

Au-Au(200 GeV) 

0.92 

0.96 

1.05 

1.06 

1.15 

1.18 


Furthermore, the results of light (anti)nuclei average 
transverse momentum {pt) calculated by PACIAE+DCPC 
model in central rapidity Pb-Pb collisions at .^/snn = 
2.76 TeV are given in Table IIIII As a comparison, the 
results [13 of light (anti)nuclei average transverse mo¬ 
mentum {pt) in Au-Au collisions at y^SNN = 200 GeV are 
also showed. Here we can see that the average transverse 
momentum of light anti-nuclei is slightly larger than that 
of the corresponding light nuclei. And the results of light 
(anti)nuclei average transverse momentum (pt) in Pb-Pb 
collisions are larger than one of Au-Au collisions. 

In Fig. m we show the integrated yields dN/dy of d, 
d, ^H, ^He and ^He divided by Apart, as a func¬ 
tion of Apart, respectively. All data points are normal¬ 
ized to the values obtained in the peripheral collisions 
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(50%-80%) with participant numbers 36. We choose 
the 50%-80% centrality bin as the reference of periph¬ 
eral collisions because of the limited statistics of IVpart 
and yield, as well as the strong fluctuation in the 80%- 
100% centrality bins. It shows that the yields of light 
(anti)nuclei and (anti)hyper-triton per participant nu¬ 
cleon increase rapidly with the increase of the number of 
participant fVpart- Obviously, this distribution properties 
of light (anti)nuclei and (anti)hyper-triton production in 
Pb-Pb collisions at .^/snn = 2.76 TeV mainly depend on 
their mass number, i.e., the greater the mass number, 
the faster the yield increases. 
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FIG. 5: Coalescence parameters Ba as a function of A^part for 
light (anti)nuclei in Pb-Pb collisions at ^snn = 2.76 TeV. 
The open symbols represent our mode results evaluated for 
different types of (anti)nuclei. The solid points are the data 
points from ALICE measurements evaluated with deuterons, 
hypertriton and ^He [sl. [3^. IstI. [3^ . 


In heavy ion collisions, the coalescence process of 
light (antilnuclei, and (anti)hypernuclei is historically 
described [isl - ll^ by the coalescence parameter B 4 . The 
differential invariant yield is related to the primordial 
yields of nucleons and is described by the equation: 




<PNa 

(PPa 


Ba{Ep 


d^Np A 

<PPp’ ’ 


( 6 ) 


where E(PN/d?p stands for the invariant yields of nu¬ 
cleons or light (anti)nuclei, and (anti)hypernuclei, and A 
are the atomic mass number, respectively. PA,Pp denote 
their momentum, with pA = App assumed. The phase- 
space information of p(p) and A(A), which is used as an 
input for the coalescence prescription, can be reproduced 
by the hydrodynamic blast-wave model. 

The coalescence parameters Ba can be evaluated by 
comparing the invariant yields of the light (anti)nuclei 
and the primordial (anti)nucleons. Fig. 5 presents the 


'^~y/B~A as a function of Npart, which remain unchanged 
from central to peripheral collisions for our model re¬ 
sults. B 2 and a/B s cal culate d ba sed on the invariant 
yields of d, d, ^He{^He), are consistent with 

the ALICE results [33,l3l,[^|^ as Npart > 100. In ad¬ 
dition, we noticed that of is smaller than 

one of ^He{^He) even though the two nuclei have the 
same mass number. This reflects the strangeness con¬ 
tent dependence of the coalescence parameter, i.e, there 
exists an additional penalty factor due to strangeness. 

In a word, we hope utilizing the study of light 
(anti)nuclei and light (anti)hyper-triton, like their yields 
and ratios, their centrality dependence, their transverse 
momentum distributions, to seek a new clue to explore 
the new production mechanism for heavier antimatter in 
relativistic heavy ion collisions. 


IV. CONCLUSION 

In this paper, we use the PACIAE+DCPC model 
to investigate the properties of light (anti)nuclei and 
(anti)hyper-triton production in Pb-Pb collisions at 
^snn = 2.76 TeV with \y\ < 0.5 acceptances. The in¬ 
tegrated yields dN/dy of light (anti)nuclei d (d), ^Ele 
as well as are calculated by the DCPC 

model for each centrality bin. The results show that 
the yields of light (anti)nuclei and (anti)hyper-triton de¬ 
crease rapidly with the increase of centrality bins. How¬ 
ever, the yield ratios of light anti-nuclei to light nu¬ 
clei and anti-hypertriton to hyper-triton are indepen¬ 
dent of centrality, which are irrelevant to the decreas¬ 
ing yields from central to peripheral collisions. In ad¬ 
dition, the transverse momentum distributions of light 
anti(nuclei) are given in the 0-10% most central Pb-Pb 
collisions. With respect to light nuclei, the transverse 
momentum that the peak of distribution corresponded 
is nearly equal to the one of corresponding anti-nuclei. 
Our model results are in agreement with the ALICE 
experimental data. We also gained that the yields of 
light (anti)matter per participant nucleon increase lin¬ 
early with the increase of the number of participants with 
Bpart. Obviously, this distribution properties of light 
(anti)nuclei and (anti)hyper-triton production mainly 
depend on their mass number. At last we calculated and 
discussed coalescence parameters Ba of light (anti)nuclei 
and (anti)hypernuclei. 
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